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LetP(u f x)=e-'^^f(^\u>0. 



The paper studies the convergence of P(u,x) to f(x) as u— ><». The results obtained 
are generalized analogs, for the interval 0<x < °° , of known properties of S. Bernstein's 
approximation polynomials in a finite interval. 



1. With a function J(t) in the closed interval 
[0,1], S. Bernstein in 1912 associated the poly- 
nomials 



B n (t) =ib(^t«(l-ty-<f(v/n), 71 = 1,2,3,- 

v-0 \V J 



(i) 



He proved that if /(0 is continuous in the closed 
interval [0,1], then B n (t)-*j(t) uniformly, as n-* °° . 
This yields a simple constructive proof of Weier- 
strass's approximation theorem. 
More generally the following theorems hold: 

Theorem A. If f(f) is bounded in [0,1] and 
continuous at every point of [a, b], where 0^a< 
6^1, then B n (t)->j(t) uniformly in [a, b]. (See 
[6], p. 66) 1 . 

Theorem B. If j(t) is bounded in [0,1] and 
continuous at a point r, then 5 w (t)— >F(r) . (See 
[3], p. 112). 

Theorem C. If /(0 satisfies a Lipschitz-H older 
condition 

1/(0-/(0 \<c\t-t'\\0<\£l, 

then |/(0~ B n {i) |<c 2 n~ x/2 ,Ci,c 2 constants (see [7], 
p. 53 ;4). B n (t) is a linear transform of the func- 
tion f(t); for the infinite interval (0,<») we define 
an analogous transform: 



v = VI 



(2) 



1 Figures in brackets indicate the literature references at the end of this 
paper. 



We shall prove corresponding theorems of approx- 
imation for this transform; u— ><» corresponds to 
n-^co in eq 1. We also sharpen theorem B to 
uniform convergence at the point r. 

Definition: A set of continuous functions P(u,x) 
is said to converge uniformly to the value S at a 
point x=£, as u->&> if P(u n ,x n )^>S, whenever 
x rt — >f and u n -^ oo , as n— > °° . An equivalent for- 
mulation is: to any e^>0 there exists a 5(e) and an 
v(d, e) so that \P(u,x)-S\<ie for \x— f|<a and 

U>7). 

2. In this section we introduce some lemmas 
for later application. 

Lemma 1. For A>0, u>0, 



W 



S ~,<\- 2 ue\ 

|p-5|£X VI ^ 



The following identity is easily verified: 

u[ 
vV 






(3) 



(4) 



it follows that 
X 2 



This proves lemma 1. 
Lemma 2. For u> 



ntv co rt/i' 

S l>^(v-uyl=ue«. 







1 t?i 



(5) 



2 M. Kac also considered the transform (2) independently, from a similiar 
point of view. 
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By Schwarz's inequality and by eq 4 



Hence 



this proves lemma 2. 
Observe that 



thus, if ^ is a positive integer 



(6) 






tf ! tf>w 



■2u S £ = 



2 2 fu-») -V=2uS" , 
ttu vi f^ivi vi^-l VI 

2u— y ^ —t=- e u , 

u\ V2tt 

by Stirling's formula. Thus the estimate (5) is 
the sharpest possible, except for a constant factor. 
3. Theorem 1. Suppose that /(x) is bounded in 
every finite interval; \ij(x)=0(x k ) for some &>0 
as z— > oo and if f(x) is continuous at a point f , then 
P{u)J) converges uniformly tof(x) at x=£. 
Consider 



e ux {P (u ;f) -jix) } =2 {f(v/u) -j(x) } (ux)' 



1] + S =S 1 +&,say 



»!' 



and assume 


that |z— f|<5. 


Let 


max |/(3)-/(0| = m(«,0 = m(5), 
for 


|Xr-£|<«, 


then m(5)^0 as 5->0. Now 


/(»/*) -/(*) =/(*/*) -/(f) +/(f ) -/(x) , 


and 


f= — x+x— f. 


In S x 




^5, 


hence, 


from 


eq 7 



(V) 



and 



|«?/tt— f[^2S 



|/(«?/u) -j(x) | ^ m (25) +m(8)S 2m (25) . 



|^|<2m(25) JZ ^=2m(2b)e" 



Next write 



Then 

Let 
Then 



$2= X) + S = $3 +£4, say. 



sup |/(«)|=?7(5), forx^S. 



ux—vyub 01 



Applying lemma 1 with \=ud, we get 



\S 3 \<2U«+8) ^ = 2U(!+5) ^r 



Finally, assuming w(z+<5)>&, 



£ 4 =o ( s 

V^>w(2:+5) 



c 

\M>M(o;+5)-fc Mi / 






*(z+S)-fc 

We apply again lemma 1, with \=u8— &>0, 
then 



Summarizing, we find 

1 , 



P(u;f(x))-f(x)=Q m(2e>) 



^->oo. 



// 



^5 2 ' (^<5-&) 2 j 
Letting u-> oo for a fixed 5, 
limsup |P(w;/(a:))-/(z)| rg 0(m(25)),u-> oo , |a;-f | ^ 5, 

from which our theorem follows. 

It can be shown easily that uniform convergence 
at each point of a closed set D implies uniform 
convergence over the set D. A similar argument 
applies to the transform (1), thus sharpening the 
theorems A and B. 

4. Theorem 2. If j(x) satisfies the Lipschitz- 
type condition 
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I hen 



WW-/fc> I <T J£=J^ 0<*<*< . , (8) <e-y±^f l*£t=*U y ^±^ j£=«L 

y, p constants, 0<p^ 1, ^ye' ux " (u^) p , , 7 

< — r ^->~^C \ V ~ U£ \^ -7=1 

\P{u]f{x))-j{;x)\^yu-pi 2 1 by lemma 2. This proves theorem 2 for p=l. 

Now from Holder's inequality, for 0<p<l 
uniformly for 0<^< oo } as w— > oo . 

We have, for ,-1 *-\ P <*> »"/(*) I = S^I? ff !/(»/«)-/(*) | 

| P(«;i) -/(*) | g «— S («*)• Jj |/(*/t») -/(*) | g (S ^) I ") {^^f ^ Vln) -J® 1 1/P |"- 

Assume that 

then 

<7 exp ((1 — p)iu;)f — 7 s ^uxe ux ) =yu~p /2 e ux - 

This completes the proof of theorem 2. 

Let p=l, and f(x)=c— x, for O^x^c, c a positive constant; f(x) =0 for x^c. Now the condition (8) 
is satisfied. Furthermore 

P{u; /) -/(c) =P(u; f) =«-- £ -^ (<;-»/«) =i e-«" S («c-») &&. 
Let [i£e]=&, then 

t^Pfa, c) > u-^e-^ 1 X! (*-») K 

v = V- 

and, from section 2 



& Lc jLfc+l / h \ M 



Thus Now 

lira infu* P(w,c)>0. p («•/)= «-*•< V} ^1 fr v / u \ 

This proves that for p= 1 the order of the estimate _, , , _ x ^ 

in theorem 2 is the sharpest possible. We do 

not know of a similar example for p<0. For _ «> (uxY n 

Bernstein's polynomials an exact result has been \ u >Pn)=e 2Lj ^j ^-i a k e 

given by M. Kac [4]. n m (-k/u)Y 

5. Suppose that /(x) is continuous in the infinite =e xu X] &* S ~ — j ^^ 

interval (0, oo). 

L- et =e _:CM S dk exp (wz exp (—k/u)) 

3=logj,0£t£l, 

/(z)=/(log l/t)=<j>(t) is continuous in O^^l. ° 

Given e>0, we can find a polynomial X) <M* =3>«(0 Clearly for «-*•<», PCwjpJ-^,,^-*) uniformly in 

so that j <M0 -2>„ (0 |<e It follows that ( °' °° } ; furtherm <>re 

\j(x) -*.(«-) ! <e, 0<z< « . /(a;) =*„(«-«) +«,(*), \e n (x) \< f , 
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and 
Here 



Now 



P(u;f)=P(uJ-p n )+P(u,p n ). 
\P(u;f-Pn)\<*, 



hence 

\P(wJ)-M\<e+\P(u;p n )-f(x)\^e + 
\P (u; Pn ) -p n (O | + \Pn (e~ x ) -f{x) | . 
Thus, the theorem: 

Theorem 8. Iif(x) is continuous in (0, oo) then 
P(u)j)-^j{x) uniformly in (0, °°). 

6. Theorem^. Ifj(x) is r-times differentiable, 
J ir) (x) = 0(x k ) as x-xn, for some K>0, and if 
j (T) (x) is continuous at a point f, then P ir) (wJ) 
converges uniformly to/ (r) (;z) at x={. 

We write l/u=h, and introduce the notation 



Aj(vh)=f(v+lh)-j(vh). 



A 2 / (vh) =AAf(vh) =j(v+2h) - 2/0+ lh) +J(vh) . 
A r f(vh)=AA T - l J(vh) 

=§(-1)* (I) jG+kh),r*0. 

Lemma 3. We have 
Differentiation gives 

i-~*¥a(f)Vc*> • 

The lemma now follows by induction. 
It is known that 



A7W 
h r 



-f T) (v), 



where 



vh<^rj<^(v-\-r)h. 



D r Q(h-J)-j"(x) 

=«-"{s +z; J. 

where we assume that |z— f|<5. Using the same 
device as in the proof of theorem 1, we get theorem 
4. For Bernstein's polynomials see G. Lorentz 
[5], and his reference to Wigert's work. 

For p=l theorem 2 and formula (6) suggest the 
following proposition: 

Theorem 5. If f(x) is bounded in every finite inter- 
val, if it is differentiable at a point f>0, and if 
f(x)=0(x k ) for some fc>0, #^o°, then 



Let 



max 



«*{P(u;/(r))-/(r)}->o,«->». 

M+h)-J(t) 



h 



-/'«■)!=*•(*, «=/»(*), 



then /z(5)-^0 as 5— >0. We may write 

/(r+*)-/(rt=V'(fl+A«cr,*), 

|.(r,*)|£M(*)far |A|<«. 



where 
Now 



P(^;/a))-/(f)=^^S 



mn-t> 



a)+ 



(5-0 



«.(«) 



where 



1** 



Utilizing formula (6) we get 



W (|i;-ttr|^«tt |»-Mf|>«« 



Using the same device as in the proof of theorem 
1, and employing lemma 2, we can complete the 
proof of theorem 3. 

The result can be generalized to higher deriva- 
tives. We restrict ourselves here to the case that 
/"(f) exists. Thus, 

/(f+A)-/(r)=¥ / (rt+|A 2 {/ / '(i-)+€(r,A)}, 
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where 

|f(T,*)|Sl(«) for \h\£t, and i?(5)-»0, 5-»0. 
Now 

p{«;y(f)}-/(f)=^s^-'{(|-f)/'(r) 



then 



+ 



I (~f)/"(f) } +< r " f ? t I (^-f) 2 ^(«), 



2 
where 



|e,(tt)|^u(8) for 



'/ 



<8. 



(9) 



It follows from formulas 4 and 6 that 



P(«;/(r))-/(0-"£/"<r)+ 






or 



«{P(«;/)-/(r)}=r ? f/"(f) + 



?S^-'( e -ttf)V«). 



2m o ^! 



We write 



fri 



V— U£\^U& \v — Wf|>U* 

= T 1 +r 2 ,say. 



From (9) and (4), 



Hence, 



Next write 



iz-.K^urc*. 



£|r.|<|,(«>r. 



r 2 = S + S =r 3 +T 4 , say, 
and note that 

i0-fJe.(«)=/(0-/(r)-(^-f)/'(f)- 

Ks- f )' y " (r) - 

Let 

sup|/(a!)|=M(f),agf 



|T 8 |<J2M(f)+f|/'(r)|+ifV" 



(f)l 



S « 2 ? -,(«*•)', 

u£-v>ui v > 



We now employ the formula (see e.g. [2], p. 200) 

|r-t|>ft* *>! V V 3 // 



It follows that 



so that 



-"'r.=o{„«p(-I|„)^ 



Finally, in view of /(a;) =0(x fc ), we have for v^>u£ 
and fc > 2 

hence 

- o C.5,- l " ! T)= o !» i «»(» f -|*)^ 

Thus, 



-«f 



T 4 =0 <w exp 



■U*-) 



(ii) 



(10) Summarizing, from (10) and (11), 

lim sup \u{P(u;M) -/(f)} -\ f/"(f)|g 5. 

But 5 is arbitrarily small, hence the theorem: 
Theorem 6. If f(x) is bounded in every finite 
interval, if it is twice differentiate at a point 
f>0, and if for some &>0 f(x)=0(x k ), x— >oo ? then 



Analogous theorems for Bernstein's polynomi- 
als were given in [1] and [8]. 

8. In the terminology of probability distri- 
bution the Bernstein polynomial corresponds to 
the binomial distribution. The distribution func- 
tion is 

T,(";)m-ty-=F n (r); 
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the linear functional B n {t ;/) is 

V j(rln)dF n (r)=ibj(!>ln) (*Wh)— • 

JO v=Q \V / 

Similarly, P(wJ) corresponds to the Poisson 
distribution; the distribution function is 



~r VI 



--6(r,x), 



and 



f*00 CO 

P(u;f)= f(rJu)dG(r,ux)=Jl< 

J0 * o=0 



(ux) 1 



here the term with the largest weight has the 

index v^ux. In the Bernoulli polynomial the 

term with the largest weight has the index v^tn. 

If instead of a function/ we consider a sequence 

S , Si, S 2 , . . . , then to the transform (1) cor- 
responds 



S( „)f(l-t)—S.,n- 

F = \v / 



which defines the generalized Euler summability, 
and to the transform (2) corresponds 

00 rpV 

e~ x S —. S t) x-> oo 

which gives Borel's summability method. 

9. To approximate a function f(x) over the 
whole real axis, we write 



/(*)' 



f(x)+f(-x) , /(s)-/(-s) _ 



: /iW+/2W, say; 



obviously 



Now 



here 



M-x)=fidt),M-x)=-M*)- 



P(u;f)=P(u;f l )+P(u;f 2 ); 



P(u-J 1 )=e- 



2&pMm> 



81 

and if we change x into— x, u into— u, we get 

P(-u;f 1 (-x))=P(u;f l (x)), 

so that our previous results are directly applicable. 
Similarly, 

P(-u;j 2 (-x)) = -P(u;f 2 (x)); 



thus for negative values of x we need only change 
u into —u, and revert to our previous results. 

10. It follows from a well-known property of 
the Beta function that 



Oi 1 ^ 1 -^^^^ ' 1 ' 1 



',n; 



hence 



f(v/u) ; so that for any Riemann integrable function 

Similarly, at first formally 

J o P(u-J)dx=T,~j(v/ U )j o e-*Vdx= 



(9) u^ {vlu) > 

the interchange of integration and summation is 
legitimate if the series ^2f(v/u) is convergent. 
Thus, the formula 

Jj P{u-J)dx=\±j (»/«) 

is valid if both sides exist. However, it is a deli- 
cate question under what conditions 



1 



J f(x)dx. 
An extensive literature deals with this question. 



w _»co U> 
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Numerical Determination of Characteristic Numbers J 

By W. E. Milne 

This paper represents a contribution to the problem of characteristic values and charac- 
teristic solutions of ordinary linear differential equations. The problem is conceived as a 
vibration problem in x and t. The partial differential equation is then approximated by a 
difference equation in both entries. The problem is now to find those frequencies in which 
a separation in x and t takes place. This is done by finding the roots of a trigonometric 
expansion of certain order. The method is applied to a number of interesting cases that 
illustrate various types of situations encountered in problems of physics and engineering. 



An improved method for obtaining the latent 
roots of a matrix has been devised by Cornelius 
Lanczos of this Bureau. He has shown further 
how this process may be used to secure the char- 
acteristic numbers (eigenwerte) belonging to a 
boundary value problem associated with an ordi- 
nary linear differential equation. The purpose of 
this note is to present a procedure for calculating 
characteristic numbers, based essentially on this 
method of Lanczos, but modified in such a way 
as to provide a simple numerical routine for the 
computation. A number of numerical examples, 
worked in full, illustrate the procedure and give 
an indication of the accuracy attained. The 
exposition is limited to the case of differential 
equations of the second order, but the method is 
capable of extension to cases of higher order. 

1. The problem. Let 



d 2 u 



du 



L(u)=P (x) %j+Pi(x) ;g+P«(*)tt. 

The problem before us is to find those characteris- 
tic values of X for which the differential system 



. du A , 

u-\-g , =0 at x=a 

u-\-G j- =0 at x=b 
ax 



(1) 



possesses nonzero solutions in the interval 
aSx^b. 

Problems of this type arise in many different 
ways in mathematical physics. For example, 
they may occur in connection with the heat 
equation 



dt 



=i<F), 



if we assume particular solutions of the form 

V=u(x)e-™, 

or in connection with the wave equation 
b 2 V 



1 The preparation of this paper was sponsored (in part) by the Office of Naval Research. 
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